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Abstract 

Massive and massless vector particles, in the framework of gen- 
eralized Weinberg- Tucker-Hammer equations, are considered in the 
Petiau-Duffin-Kemmer formalism. We obtain solutions of equations 
in the form of matrix-dyads. The Lagrangian is defined in the matrix 
form and the conserved electric current and the energy-momentum 
tensor are obtained. The canonical quantization is performed and the 
propagator of massive fields is found in the formalism considered. 

1 Introduction 

The vector fields play very important role in the Standard Model which 
unifies weak and electromagnetic interactions of elementary particles. After 
the spontaneous breaking symmetry massless vector fields acquire masses due 
to the Higgs mechanism. But the existence of scalar Higgs bosons in nature 
is questionable and is verified at the Large Hadronic Collider. Therefore, the 
alternative description of massive and massless vector particles is of great 
theoretical interest. The second-order Proca equations, that are used in 
the theory of vector particles, can be represented in the form of the first- 
order Petiau-Duffin-Kemmer (PDK) relativistic wave equation [JJ, [2], [3] 
(see also [J]). The PDK form of equations is very convenient for different 
applications [5j. We have used the PDK form of equations in the Dirac- 
Kahler theory j6], in the field theory of multi-spin particles [7], [8], and in 
generalized electrodynamics [9]. 

In this paper, we consider generalized Weinberg- Tucker-Hammer equa- 
tions (WTH) for massive and massless vector fields, introduced by Dvoeglazov 
[10], in the PDK formalism. We obtain solutions in the form of projection 
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matrix-dyads, find the conserved electric current and the energy-momentum 
tensor, and perform the canonical quantization in the PDK form. 

The paper is organized as follows. In Sec. 2, we represent generalized WTH 
equations for massive and massless vector fields in the form of generalized 
PDK equations. Solutions of the wave equation for a free fields are obtained 
in the form of matrix-dyads in Sec. 3. We define the Lagrangian in the matrix 
form and find the conserved electric current and the energy-momentum tensor 
in Sec. 4. The canonical quantization of fields is performed and the propagator 
is obtained in the PDK form in Sec. 5. We discuss the results obtained in 
Sec. 6. Some useful matrix relations are found in Appendix. 

The Euclidean metric is used and the system of units H = c = 1 is 
explored. Greek and Latin letters run 1,2,3,4 and 1,2,3, respectively, and, 
we imply the summation in the repeated indexes. 

2 Generalized WTH equations in the PDK 
formalism 

Let us consider generalized WTH equations for vector fields, introduced by 
Dvoeglazov [10J: 

(lapd a d p + Ad 2 a - Bm 2 ) V>(z) = 0, (1) 

where, 7 Q/ g are 6 x 6-matrices given in [TT], = (d/dx m ,d/d(it)), A, B are 
dimensionless parameters, and m is the mass parameter. The wave function 
ip =column(E+zB, E— zB), (where E a = iF a 4, B a = (l/2)e a bcFb c ) realizes the 
(1, 0)©(0, 1) representation of the Lorentz group. Eq.(l) is the generalization 
of the Weinberg [T2] and Tucker-Hammer equations [13J. At A = 0, B = 1 
one arrives at the Weinberg equation, and at A — 1, B — 2, we have the 
Tucker-Hammer equation [10J. In the component form Eq.(l) is equivalent 
to the second-order equation for the antisymmetric tensor F a p [TO] : 

d^Ftf - dpdpF^ + -^8%F a p - ^m 2 F a/3 = 0, (2) 

We study Eq.(2) in the PDK form as for the massive case, m ^ 0, as well as 
for the massless case, m — 0. 
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2.1 Massive vector fields 

We introduce the 10-component wave function 



= {^(x)} 



(A = |i,H)- (3) 



With the help of the elements of the entire matrix algebra e A ' B , with products 
and matrix elements 



/ 'M,N\ r r „M,A^tSA\l 

where A, B, M, N — /i, [/iv], Eq.(2) can be represented in the form 



(4) 



M _|_ gM." 



me 



+ 



d. 2 , + 

2,771 2 



(5) 



*b(x) =0. 



We imply a summation over all repeated indices. Defining the 10 x 10 matrices 

(6) 



^ _ £ f,M _|_ £ M,^ p — £ M>^ p — I £ M,M ; 



Eq.(5) takes the PDK matrix form 

I- A 



Ppd,* + mP + P 



B 



9J, + —m 
2m * 2 



= 0, 



(7) 



where the matrices /3 M are Hermitian matrices, /?+ = (3^. Eq.(7) is the gen- 
eralized PDK equation. The projection operator P = P + extracts the four- 
dimensional vector subspace (ip^) of the wave function and the projection 
operator P = P + extracts the six-dimensional tensor subspace corresponding 
to the VVi/] - The matrices /3 M obey the PDK algebra 



and matrices P, P are projection matrices 

p 2 = p ) p2 = p p + p = x pp = pp = 0? 



(8) 



(9) 
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The Lorentz group generators in the 10-dimension representation space 
are given by [7J, [5] 



(10) 



and obey the commutation relations 

\Jpai J pv 



^apJpv ^pyJcrfi ^ ppJ av ^auJppi 



[A\> JpA = §\pP>v - hvP, 



I-'- 



p j 

1 i ° pv 



0, [P,J, 



pu\ 



0. 



(11) 



Eq.(7) is form-invariant because of Eq.(ll). The Lorentz-invariant is = 
fy + 7]ty, where is the Hermitian-conjugated wave function, and the Her- 
mitianizing matrix rj is 



n 



2 



(12) 



The f] is the Hermitian matrix, 77 + = rj, and obeys the relations: rj/3 m = 
—f3 m f] (m=l,2,3), 77/5 4 = /3 4 r]. With the help of these relations, one finds the 
"conjugated" equation 



Pfir -mP- (U 1 -^ + |m ) P 



0. 



(13) 



where \I/ = (ip*, — V'jwi) an d the complex conjugation * does not act on the 
metric imaginary unit i of fourth components of the wave function, ip* = 

(ip^, i'lpo) , and so on. 



2.2 Massless vector fields 

To get the massless WTH equation, we put m = in Eq.(2) and arrive at 



d a duFpB — dRdnFu a H — d u F a s — 0. 



(14) 



In the same manner as for the massive case, we introduce the 10-component 
wave function 



K ^u-Ppu (^) 
■Ppvi^j 



(15) 



where k is the new parameter with the dimension of the mass. We need the 
dimensional parameter to have the components of the wave function 
with the same dimension. Physical values do not depend on this parameter, 
and k is absorbed by the renormalization of fields [9]. As a result, Eq.(14) 
takes the form similar to Eq.(7) for massive fields 



13,8, + kP + P-dl 

K 



V{x) = 0, (16) 



where C is an arbitrary parameter (C = (A — l)/2). Thus, the equation 
of motion and other equations for the massless fields can be obtained from 
the massive case by the replacements m — > k, B = 0. Therefore, in the 
further consideration, we write down only the expressions for the massive case 
implying that they are valid also for the massless fields after the replacements 
m — > k, B = 0. Eq.(14), (16) represent the modified Maxwell equations. 

3 Solutions of the matrix equation 

In the momentum space, for the positive (+p) and negative (— p) energies, 
Eq.(7) reads 

(±ip + mP + \P)m(±p) = 0, (17) 
where p = P^p^, the four-momentum is p M = (p, ip ), and 

A- 1 9 B 

p 2 = p 2 — Pq. The matrix of Eq.(17) 

A± = ±ip + mP + \P, (19) 

obeys the minimal polynomial (see Eq.(A3) in Appendix): 

(A± - m) (A± - A) [(A± - m) (A± - A) + p 2 ] = 0. (20) 

There exist non-trivial solutions of Eq.(17) if det A± = or the eigenvalue 
of the A± is equal to zero. If A ^ 0, this requirement leads to the dispersion 
relation 

p 2 + Am = 0. (21) 
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It follows from Eq.(20) that there are also other eigenvalues of the matrix 
A±: m and A. Taking into account Eq.(18), one obtains from Eq.(21) 



9 m 2 B , s 

p + aTT = (22) 



Thus, the mass of the vector field is 



M = m Vlfr (23) 

Another case which leads to det A± = 0, is A = 0. Then Eq.(18) results at 
A = to the physical mass 



M ' = m VI^T (24) 

The formulas (23), (24) are in the agreement with results obtained by Dvoeglazov 
[TO] . From Eq.(23), we find the restriction on the parameters A and B: 
B/(A + 1) > 0, and in the case of Eq.(24), B/(A - 1) > 0. It follows from 
Eq.(23) that for Weinberg [TJ (A = 0, B = 1) and Tucker-Hammer [JJ 
(A = 1, B = 2) equations, we have the same physical mass of vector par- 
ticles, M = m. The case (24) is not realized for Tucker-Hammer equations 
[T3] [A — 1, B — 2, A ^ 0). Therefore, we imply that A ^ and the physical 
masses of vector particles are given by Eq.(23). On-shell, when Eq.(22) is 
valid, and m^A the minimal polynomial equation (20) reduces to 

A± (A± - A - m) (A± - m) (A ± - A) = 0. (25) 

Eq.(25) allows us to obtain solutions of Eq.(17), A^{p) = 0, in the form of 
the projection matrix [14]. Indeed, the matrix 

n± = (A± - A - m) (A ± - m) (A± - A) 

• . \ ( 26 ) 

= ^fNip l±impP ± iApP — Am] , 

where A^ is the normalization constant, obeys the equation A-i-H-i- = 0. It 
means that every column of the matrix n_|- is the solution to Eq.(17). The 
requirement that the n_|- is the projection matrix 

ni = n± , (27) 
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gives the value 

N = ~^V V- (28) 

Xm(X + m) v ' 

Eq.(28) can be verified with the help of the minimal Eq.(25). It should be 
noted that Eq.(25) represents the minimal polynomial for the non-degenerate 
case m ^ X (see Eq.(A3) in Appendix). But the case m = X is realized as 
for the Weinberg equation (A = 0, B = 1) as well as for the Tucker-Hammer 
equation (A = 1, B = 2). This follows from Eq.(18),(23). For this special 
case, m — X, one should explore the minimal polynomial equation on-shell 
(see Eq.(A4) in Appendix) 

A± (A ± - m) (A± - 2m) = 0. (29) 

which corresponds to the PDK equation for spin-1. Indeed, if m = X, Eq.(17) 
becomes 

(±ip + m)V(±p) = 0, (30) 

as P + P = 1. But Eq.(30) is the PDK equation for vector fields and the 
solution in the form of the matrix-dyad to this equation is well known [14]. 
We find the projection operator from Eq.(29) 

T PDK _ 1 / A 2 o^A ±?(±»P-m) 



II— = _ - 3mA ± + Wj = (3D 

Eq.(26) converts to Eq.(31) at m = X. 

The projection operators extracting spin projections ±1 and are given 

by m, x 

S(±i) = 2°p(°p ± !) > S( Q ) = l-o%. (32) 

where the spin operator is 

°V = _ 771 i^abcPaJbc = ~] \tabcVaPbPc (33) 

2 I P I I P I 

and commutes with the matrix of Eq.(17) A±: [A_|-,crp] = 0. The projection 
operators satisfy the relations: S?±i) = ^(±1); <3(±i)S(o) = 0, = S(p). 
The operators (32) commute with the mass projection operator (26) and, as 
a result, we obtain projection operators extracting solutions to Eq.(17) for 
spin projections ±1, (for states of particles with the mass m) in the form 
of matrix-dyads 

n±%i) = *±i ■ *±i, n±s {0) = • *o- (34) 
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The matrix elements of the matrix-dyad is • ^)ab — ^a^b- For the 
case of the Weinberg and Tucker-Hammer equations, one has to make the 
replacement in Eq.(34) IT ± -»■ T1^ DK . Eq.(34) can be used for different elec- 
trodynamics calculations of processes with vector particles in the framework 
of generalized WTH model. 

For the case of massless WTH equation, the parameter on-shell is A = 
(p 2 = 0). Then the operator of the equation becomes 

A± = ip + K P, (35) 

and obeys the minimal matrix equation on-shell (see Eq.(A6) in Appendix) 
as follows: 

A± (A± - k) 2 = 0. (36) 

As pointed in [15], in this case there are difficulties to obtain the solutions in 
the form of projection operators because the zero eigenvalue of the operator 
A± is degenerated. In the case of Maxwell equations such difficulty was 
overcome in [9] by using the general gauge and adding a scalar field. 



4 Lagrangian and conserved currents 

The Lagrangian of generalized WTH model in the PDK formalism is given 
by 

2 v ; v 2 J (37) 

One may construct another Lagrangian by adding the four-divergence to 
Eq.(37) that does not change equations of motion. It can be verified that 
Euler-Lagrange equations following from Eq. (37) give the equations of motion 
(7), (13). We obtain also the Lagrangian (37) in terms of fields if) a 



(38) 
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mB I - A 

+— ^W]^W] 



+ c.c. 



2 ^ym^ym 2m 

where the c.c. is the complex conjugated expression. The electric current 
density is given by 



dC 



dC 



x. 



(39) 



Replacing Eq.(37) into Eq.(39), we obtain the electric current density 

I — A 



M(x)pj!(x) + \^(x)Pd li 9(x) - (d»V(x)) PV(xj\ . (40) 

2m L v ' 



One may check with the help of equations of motion (7), (13) that the electric 
current is conserved, d^j^x) = 0. The second term in Eq.(40) vanishes 
at A = 1, corresponding to the Tucker-Hammer equation, and j^{x) takes 
the form similar to the Dirac theory. Using the general expression for the 
canonical energy-momentum tensor 



T, 



dC 



fW 



d(d^(x)) 



dvV(x) + d v V(x) 



dC 



d(d^(x] 



(41) 



we obtain from the Lagrangian (37) 



T, 



fJbV 



~ (d^(x)) 0^(x) - U>{x)(3,d^{x) 



+ 



I- A 



2m 



2 



PV(x) + ^(x)Pd tl d v ^{x) 



(42) 



One may verify that the canonical energy-momentum tensor is conserved 
tensor, d^T^ u = 0. At A = 1 Eq.(37),(42) are simplified. It should be noted 
that the Lagrangian C, Eq.(37), does not vanish if A ^ 1 for fields satisfying 
the equations of motion. It is easy to write down expressions for and T MJ/ 
in the component form exploring Eq.(3),(4),(6). 



5 The canonical quantization 

Solutions to Eq.(7) with definite energy and momentum in the form of plane 
waves are given by 

*i ±} (*) = t/^^*.(±p) exp(± W ), (43) 
V 2 PoV 
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where V is the normalization volume, p 2 = p 2 — p^ = — m 2 , and s is the spin 
index (s = ±1,0). The 10-dimensional function \I/ S (±p) obeys Eq.(17). We 
use the normalization conditions 



v 



^f\x)P^f\x)d z x = ±8 ss ,, [ V i s ± \x)fa¥j\x)d 3 x = 0, (44) 



where ^^(x) = (y^(xj^ rj. Eq.(44) correspond to the normalization on 
the charge. In the second quantized theory, the field operators are represented 

^{x) = Y J [a P ^ { +\x) + b + p ^\x)_ , 



p,s 



p,s 



(45) 



where the positive and negative parts of the wave function are given by 
Eq.(43). The creation and annihilation operators of particles, a p s , a PtS , and 
antiparticles, b ps1 b PtS , obey the commutation relations as follows: 

[a PjS , ap jS ,] = 5 ssl 8pp>, [a PyS ,a p/)S >] = [a+ s ,a+, jS ,] = 0, 

[bp,; K',s>\ = $ss'Spp>, [bp, s , bp>, s '} = [6+ s , 6j jS ,] = 0, (46) 
K, s , b p , :S ,\ = [a p>8 , 6j >a ,] = [a+ s , 6 pV ] = [a PjS , 6j >a ,] = 0. 
We obtain the commutation relations for different times from Eq.(43)-(46) 

V M (x),V N (x?)] = \y M (x),y N (x')] = 0, [v M (x),V N (x')] = N MN (x,x'), 

N MN (x,x') = N^ N (x,x') - Nm N (x,x'), 

N^ n (x,x>) = J2(¥ s + \x)) m (¥;\x>) 

p,s ^ 

One finds from Eq.(43),(47) 

N£n(*, = E ^ (^(±p))m {*s(±p)) n exp[±ip(x - *')]• (48) 



(47) 



iV 



iV 
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It follows from Eq.(32) that the equation 

S i+ i) + 5 ( _d + 5 ( o) = 1, (49) 
holds, and result, we obtain from Eq.(34) 



, . ±ip ( ±impP ± iXpP — Am 

? • (».<*>)) = n ± = 1 Am(m + A) 

Taking into account Eq.(50), one finds from Eq.(48): 



(50) 



If ±ip {±impP ± iXpP — Am . . 
N MN (x, x') = ^7 I " ^ ] exp[±ip(x - x')} 



' ' Pudn (m^Pd/j, + XPf.Pd^ - Am 
Am 

where the singular functions are given by [16] 



MN 

(51) 



MN 



A+(x — x') 



p *PoV p ^Pqv 

It follows from Eq.(18) that the A (at A ^ 1) is the operator because in the 
configuration space p 2 = —d 2 . With the help of the function [16] 

A (x)=*(A + (x)-A_(x)), (53) 

from Eq.(47),(51),(53), we obtain 

iV M 7v(x,x') = -^ M^ \m Ao(x-x'). (54) 

V / M AT 

It should be noted that the function Aq(x) vanishes when x 2 = x 2 — t 2 > 
[T6] . Using Eq.(8), one finds the relation 

/ _ \ Ppdp (m^Pdf, + XppPdu - Am) 
+ mP + XP) Am 

(55) 
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(mP + XP) 0,d, 

— ( — Am 



Am 



Taking into account the equation [16] (d 2 — Xrnj A±(x) = 0, and Eq.(51), 
from Eq.(55), we arrive at 

(^ + m? + AP)iV ± (i,i') = 0. (56) 



For bosonic fields the chronological product of two operators is defined as 

{ ^_ M (x 1 )^ N {x 2 ), h > t 2 , 1 

I ^n(x2)^m(xi), t 2 >t 1 | ' 



T^ M (x 1 )^ N (x 2 ) 



Then the vacuum expectation of chronological pairing of operators (the prop- 
agator) becomes 

(T^ M (x)^ N (y)} = N c MN (x - y) 



= 0(x - y ) N^ N (x -y) + 9(y - %o) N MN (x - y), 
where 6(x) being the theta-function. We obtain from Eq.(57): 



(57) 



_ (3 u d u (m3 u Pd u + Xf3 u Pd u - Am j 

(T9(x) ■ *(y)> = " " { \ - U c (x - y), (58) 

Am 

where 

A c (x -y) = 9(x - y Q ) A + (x -y) + 9(y - x ) A_(x - y). (59) 
Using the equation [16] {d 2 , — Xrnj A c (x) = i8(x), one finds from Eq.(55),(58) 



^ + mP + XP) (T*(x) ■ V(y)) = % W(x - y). (60) 



mP + XP 
Am 



At A 7^ 1 (m 7^ A), we have non-local operators in Eq.(58),(60) because the 
variable p 2 = —d 2 enters the A. Expressions for PDK theory follow from 
Eq.(58),(60) at m — A (A — 1, B — 2), and as a result, there are only local 
operators. 
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6 Discussion 



We have considered the generalized WTH equations, that describe vector 
fields in the PDK formalism, Eq.(7). Eq.(7) is valid for any parameters A and 
B. For the Tucker-Hammer equation at A — 1, B — 2 Eq.(7) reduces to PDK 
equation. For the Weinberg equation, at A = 0, B = 1, we obtain from Eq.(7) 
the generalized PDK equation which is more complicated because it contains 
the second derivatives. To consider solutions to the wave equation (7) for 
particles in external electromagnetic fields, one has to make the replacement 
dp — > d^ — ieA^ (A^ is the vector-potential of an electromagnetic field). Then 
solutions to the Weinberg equation will be different compared to solutions 
of the PDK equation. Therefore, physical observable also will be different. 
The equation for the massless particles, Eq.(16), is also different (at C ^ 0, 
A ^ 1) comparing with the Maxwell equations in the matrix form [15], [9]. 

Solutions of equations obtained in the form of matrix-dyads allow us to 
make calculations of processes in the quantum theory of vector particles, 
described by generalized WTH equations, in the simple manner. From the 
generalized PDK Lagrangian, we have found the conserved electric current 
and energy-momentum tensor in the matrix PDK form that are convenient 
for analyzing. Quantization of fields and the propagator obtained make it 
possible to use the perturbation theory for different quantum calculations. 

7 Appendix 

From Eq.(9),(19), we obtain 

A| = -p 2 + (m + A) A± - Am, (Al) 

A| - (to + A) A| + AmA± = ^ip 2 p - mp 2 P - Xp 2 P, (A2) 
where 

fp = pf = p 2 £ ^ - Vil p v e v ^, fP = Pf = p^p u e M > [x " ] . 

Squaring Eq.(Al) and using the relation p 3 = p 2 p, after some algebraic ma- 
nipulations, we find the minimal polynomial of the matrix A± 

(A± - m) (A± - A) f(A± - m) (A± - A) + p 2 ] = 0. (A3) 
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On-shell, one has to put p 2 = —Am in Eq.(A3). Eq.(A3) holds for any m and 
A. At a particular case, at m = A the equation is simplified and becomes 



(A± — m) (A± — m) 2 + p 2 



0. 



(AA) 



It should be noted that if we put m = A in Eq.(A3), we get Eq.(A4) multiplied 
by the factor (A± — m), but the real minimal equation is (A4). For the 
massless case the matrix (35) obeys the equation as follows: 



A± (A± - k) [A± (A± - k) + p 2 
and on-shell p 2 = 0. 



0. 



(A5) 
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